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In high-temperature superconductivity, the process that leads to the forma- 
tion of Cooper pairs, the fundamental charge carriers in any superconductor, 
remains mysterious. We use a femtosecond laser pump pulse to perturb super- 
conducting Bi2Sr2CaCu20 8+ 5, and study subsequent dynamics using time- 
and angle-resolved photoemission and infrared reflectivity probes. Gap and 
quasiparticle population dynamics reveal marked dependencies on both exci- 
tation density and crystal momentum. Close to the J- wave nodes, the super- 
conducting gap is sensitive to the pump intensity and Cooper pairs recombine 
slowly. Far from the nodes pumping affects the gap only weakly and recombi- 
nation processes are faster. These results demonstrate a new window into the 
dynamical processes that govern quasiparticle recombination and gap forma- 
tion in cuprates. 

The lifetime of Bogoliubov quasiparticles, the low energy excitations of a superconduc- 
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tor, contains a wealth of information pertinent to the origin of superconductivity in a given 
material (1). This lifetime reflects two distinct processes: quasiparticle scattering and recom- 
bination. In the former, a quasiparticle scatters from one momentum state to another, conserv- 
ing the fermionic particle number. Recombination, on the other hand, refers to interactions in 
which two quasiparticles annihilate. To conserve energy and momentum, recombination must 
involve emission of other excitations, for example phonons or magnons, to which the quasipar- 
ticles are strongly coupled. Measurement of quasiparticle recombination rates as a function of 
their energy and momentum can, in principle, provide direct information about the interactions 
that induce Cooper pairing and superconductivity (2). Only very recently, with the demonstra- 
tion that angle-resolved photoemission spectroscopy (ARPES) can be performed with ultrashort 
laser pulse sources (3-9), have measurements with the necessary energy, momentum, and time 
resolution become possible. We present the results of experiments that use synchronized laser 
pulses to perform time-resolved ARPES measurements of quasiparticle recombination and gap 
dynamics in the high-temperature superconductor Bi2Sr2CaCu20 8+ §. 

Measurements were performed at 18 K on an optimally doped sample with a critical temper- 
ature (T c ) of 91 K. A transient state is created with an infrared laser pump pulse (hv = 1.48 eV) 
and measured via photoemission shortly thereafter, with a temporal resolution of 300 fs, using 
an ultraviolet probe pulse (hv = 5.9 eV). The experiment benefits from high momentum and 
energy resolution (0.003 A -1 and 23 meV, respectively) and the ability to explore low pump 
fluences (2-15 jUJ/cm 2 ). 

Figure 1 shows typical equilibrium and transient ARPES dispersions (t = — 1 ps and t = 0.6 
ps respectively) for cuts along nodal and off-nodal directions in fc-space. Here the time origin 
t = coincides with the application of the pump pulse. The off-nodal cut has an equilibrium 
gap of 15 meV. In both cuts a well-defined kink (marked by arrows in panels (A-B) and (F-G)) 
(10, 11) separates sharply defined coherent dispersive features from poorly defined incoherent 
features, as also visible in the selected energy distribution curves (EDCs) shown in panels (D) 
and (I). The following changes are evident in the transient spectra: 1) a decrease of intensity 
below the Fermi level (Ep) and slight broadening in the coherent spectra (panels (C-E), and 
(H-J)), similar to a previous report for nodal quasiparticles (5) and mainly confined below the 
kink binding energies (10, 11); 2) an overall transfer of spectral weight across Ep (panels (C), 
(E), (H), and (J)), indicating the creation of transient quasiparticles; and 3) a small shift of the 
spectral peak toward Ep in the off-nodal cut (panels (H) and (I)), indicating a partial closure of 
the superconducting gap. 

Figure 2 shows the temporal evolution of the superconducting gap in response to photoex- 
citation, as extracted from symmetrized EDCs at kp, the Fermi wave vector (12). Panels (A-B) 
and (C-D) correspond to two representative cuts at = 32° and § = 27°, respectively, with 
defined according to the inset of panel (E). These spectra indicate very different responses of the 
gap amplitude to photoexcitation for the two cuts. The gap is relatively insensitive to fluences 
below 5 /iJ/cm 2 but 13 /iJ/cm 2 induces a clear reduction in size. As shown in panel (E) the 
gap closer to the node decreases by 55% of its equilibrium magnitude, while the gap at § = 27° 
decreases by only 20%. This may indicate different dynamics inside and outside the Fermi arc, 
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Figure 1: Typical ARPES dispersions before and after pumping for nodal (0 = 45°) and gapped 
(0 = 31°) regions of £-space. The incident pump fluence was 5 /iJ/cm 2 . (A) Equilibrium 
(t = — 1 ps) and (B) transient (t = 0.6 ps) energy-momentum maps for the nodal state. Data are 
shown with identical color scales. The inset shows the location of the cut. The arrow marks the 
position of the dispersion kink. (C) Subtraction between (A) and (B). Blue indicates intensity 
gain and red intensity loss. (D) Energy distribution curves (EDCs) from k\ to for equilibrium 
(in black) and transient (in red) states. EDCs are shifted vertically for ease of comparison. (E) 
Difference between transient and equilibrium EDCs, integrated across the double black arrow 
in panel (C). (F-J) Same as (A-E) but for a gapped (off-nodal) momentum cut. Spectra have 
been corrected for detector non-linearity. The diagonal line in the lower right portion of (F) and 
(G) is the edge of the detector. 



3 



< 




-0.04 0.04 -0.04 0.04 -0.05 0.05 -0.05 0.05 
Energy (eV) Energy (eV) Energy (eV) Energy (eV) 



1.0 
0.9 



— 0.8 



< 0.7 
0.6 
0.5 
0.4 



(|) = 27 D , 15 nJ/cm 
<|) = 32 . 13 nJ/cm 2 



J_ 



J_ 



.a 
< 



F i» 



• <|> = 27°, 15nJ/cm 
<|> = 32°, 13 nJ/cm 2 



_L 



I I L 



2 4 

Delay (ps) 



2 4 6 
Delay (ps) 



t = 3 ps 



t = -0.9 ps 



Figure 2: Evolution of the superconducting gap following pump excitation. Symmetrized EDCs 
at kp for § — 32° at low (A) and higher (B) fluence. The gap is obtained by fitting to a phe- 
nomenological model (72), but can be approximated by halving the distance between positive 
and negative peaks. Bold curves correspond to t = 0. For additional gap fitting details, see 
supplementary information online (15). (C-D) Analogous EDCs for a cut at <j) = 27°. (E) Gap 
magnitude normalized by its equilibrium value vs. pump-probe delay for momentum cuts at 
= 27° and = 32°. (F) Gap magnitude, inverted and normalized by maximal change upon 
pumping in order to compare recovery rates. 
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which is reported to end rather abruptly at = 30° for samples of this doping (13, 14), although 
studies farther from the node are needed. Gap recovery rates are illustrated in Fig. IF, where the 
curves from panel (E) have been inverted and rescaled by their maximum change. The initial 
recovery rate is slower for states closer to the node (0.9 ± 0.6 ps _1 ) than for states farther from 
the node (1.3 ±0.6 ps 1 ), although the contrast is less apparent than that between amplitudes. 

Figs. 3 and 4 show quasiparticle recombination dynamics. In this low fluence regime the 
gap is almost unchanged for most of the recovery process, so quasiparticle recombination is 
largely decoupled from gap dynamics. Fig. 3A-C shows the temporal evolution of the above- 
Ep spectral change AI for representative nodal and off-nodal fc-space cuts, where AI is defined 
by the integrated intensity change across the blue and black double arrows in Fig. 1C. The 
spectral change is nearly symmetric above and below £> in this fluence regime, so we focus on 
the intensity above Ep because of its superior statistics and smaller background. Faster decay 
rates occur at higher fluences and off-nodal momenta, an effect which cannot be explained by 
equilibrium heating (15). 

Figure 3D summarizes the dependence of quasiparticle recombination on fluence and mo- 
mentum. The rate fo is defined by fitting the decay curves at short times to the convolution of 
a Gaussian and decaying exponential (16, 17). In line with Fig. 3A-C, two prominent decay 
rate trends are apparent: 1) fluence dependence, with faster initial decay rates yo occurring at 
higher fluences; and 2) momentum dependence, with off-nodal decay rates increasing faster 
with fluence than those at the node. The first trend implies that intrinsic quasiparticle recom- 
bination processes are observed (17, 18). The second trend indicates that this recombination 
occurs more rapidly in off-nodal regions of k- space than at the node. The fluence dependence is 
also complementary to ultrafast studies using all-optical techniques, which report a dramatic de- 
cay rate fluence dependence, particularly in the low fluence regime (16-20), and there is overall 
agreement between the ARPES results and a time -resolved reflectivity measurement taken on 
the same sample (gray circles in Fig. 3D). The decay rate measured by reflectivity is uniformly 
faster than nodal ARPES decay rates (cuts 5-6), but slower than the off-nodal rates (cuts 1-3), 
suggesting that optical spectroscopy provides an effective momentum-integrated average of the 
quasiparticle population. Interestingly, along both § = 27° and § = 32° directions quasiparti- 
cle decay rates are slower than the gap recovery rates in Fig. 2: compare 1.3 ±0.4 ps -1 and 
0.9±0.6ps 1 for the gap recovery vs. 0.49 ±0.06 ps -1 and 0.37 ±0.01 ps -1 for the intensity 
recovery. This indicates that the superconducting gap recovers well before the nonequilibrium 
quasiparticle population drops to zero. A potentially related effect occurs at equilibrium in the 
BCS model, where the gap becomes large for T only slightly below T c . Finally, we note that 
the fluence and momentum dependencies reported here are in contrast to the findings of an- 
other time-resolved ARPES work (6), where the quasiparticle recombination rate was reported 
to be independent of both fluence and momentum. The discrepancy might be explained by the 
higher fluences used in the previous work, which likely result in a complete closure of the su- 
perconducting gap, or by the coarser momentum and energy resolution compared to the present 
study. 

As noted above (see also supplementary information online (15)), the fluence dependence 
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Figure 3: Quasiparticle recombination dynamics vs. pump fiuence and crystal momentum. 
ARPES data correspond to intensity change above Ep (A/) as integrated between the blue 
and black double arrows in Fig. 1C. Time -resolved reflectivity rates correspond to fractional 
change in reflectivity. (A) Nodal decay curves at 2, 5, and 10 /iJ/cm 2 , normalized to the same 
amplitude. (B) Analogous off-nodal decay curves (0 = 31°). (C) Overlay of nodal and off- 
nodal curves at the same fiuence. (D) Initial decay rate yo vs. fiuence, obtained by fitting decay 
curves at short times (for AI(t) > A/o/2) to the convolution of a Gaussian and the function 
f(t) = A/oe~ %< - f ~ ? °)®(/ — to), where AIq and to are additional fit parameters. Time-resolved re- 
flectivity rates were multiplied by 3/2 in order to take the finite penetration depth of the optical 
probe into account (15). 
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of 70 means that ARPES decay rates are connected to intrinsic quasiparticle recombination 
processes. Time-resolved optical measurements indicate that the total (momentum-integrated) 
population of photoexcited quasiparticles n ex {t) can be described by a bimolecular rate equation 



where R is a quasiparticle recombination constant, and nj is the population of thermal quasipar- 
ticles. This is a special case of the Rothwarf-Taylor model of quasiparticle recombination (21), 
which has been successfully used to model dynamics of both conventional and high-temperature 
superconductors (17-20, 22-28). The general model also incorporates negative feedback from 
"hot" bosons (hco > 2A) that are created by the quasiparticle decay process. For low fluence 
and temperature the model reduces to Eq. 1 under both weak feedback and strong feedback 
(boson bottleneck) scenarios, which beget differing interpretations of the coefficient R (17, 23). 
In both approximations bimolecular recombination is the active ingredient in fluence-dependent 
dynamics. 

In contrast to time-resolved optics, ARPES measures the momentum-dependent nonequi- 
librium quasiparticle density n k (t). The short-time fluence-dependent recombination dynamics 
are given by 



where R kk / is a modified recombination coefficient for the interaction between quasiparticles 
at specific points k and k 1 in reciprocal space. A weighted average of R kk t over k' is given by 
dyo/dF, the rate of increase of the initial decay rate yo with fluence (15). Fig. 4 shows an 
analysis of dyo/dF as calculated by fitting straight lines to the data in Fig. 3D for fluences 
F < 12 /iJ/cm 2 . It is clear that dyo/dF increases with decreasing Fermi surface angle 0, which 
means that the rate of recombination is enhanced as the quasiparticle momentum moves farther 
from the node. 

One potential scenario for the momentum dependence of the recombination rates is that with 
increasing distance from the node the quasiparticle energy and momentum approach resonance 
with charge or spin density wave fluctuations to which the electrons are strongly coupled. For 
example, a prominent neutron spin resonance is observed in Bi2Sr2CaCu20 8+ g along the (1,1) 
momentum vector (29). Resonance between this mode and a quasiparticle pair would occur 
at a Fermi surface angle of about 12°, leading to the prediction of a peak in dyo/dF at this 
Fermi surface angle. We believe that demonstrating that recombination can be mapped using 
time-resolved ARPES and observing its strong momentum dependence will further stimulate 
development of pulsed sources that are capable of reaching all the relevant regions of momen- 
tum space. 
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Figure 4: Initial rate of increase djo/dF as extracted from straight line fits to the data in Fig. 
3D for fluence F < 12 /iJ/cm 2 . The horizontal axis corresponds to the Fermi surface angle. 
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Materials and Methods 

The time-resolved ARPES experiment uses the same setup as that reported in Ref. (5). We 
induce a nonequilibrium state by exciting the sample with an infrared pump pulse (hv = 1.48 
eV). This photon energy is sufficient to energize the system's quasiparticles, but does not eject 
them from the sample. Shortly after photoexcitation, the sample is excited with a probe pulse 
in the ultraviolet (frequency-quadrupled from the pump using two BBO crystals, resulting in 
hv = 5.9 eV). The second pulse photoemits the electrons, and their energies and momenta are 
subsequently measured using a hemispherical analyzer with a 2D imaging detector (SPECS 
Phoibos 150). The total experimental energy resolution, including probe beam bandwidth, was 
measured to be 23 meV. The momentum resolution at a given photoelectron kinetic energies 
and emission angle was 0.003 A -1 . Time resolution was achieved by a motorized translation 
stage which adjusts the relative pump and probe path lengths, and therefore the time delay 
between the two. The cross-correlation between the two pulses, measured using the response of 
polycrystalline gold 1 eV above the Fermi level (Ep), was about 300 fs. The laser repetition rate 
was 543 KHz. The typical pump pulse beam profile was 100 /im at full-width half-maximum 
(FWHM), and the typical probe pulse profile was 40 /im FWHM. The sample was mounted 
on a 6-axis manipulator, cleaved and maintained at 18 K at a pressure below 5e-ll Torr, and 
oriented by taking ARPES Fermi surface maps. ARPES spectra were corrected for detector non- 
linearity and obtained in a photocurrent regime where space charge effects were small. Spectra 
at different delay times were also normalized to each other at high binding energy (100-200 
meV), where the effect of photoexcitation is not detectable (5), but where the overall intensity 
is seen to vary slightly with time due to variations in photocurrent. The typical correction, over 
the course of 20 ps, was by less than 2% of the high binding energy intensity. 

The time-resolved reflectivity experiment was performed using the same laser system as 
the ARPES experiment, also operating at a repetition rate of 543 kHz. Both pump and probe 
pulses were infrared and were focused onto a 54 /im spot on the sample. The pump delay 
was modulated using a rapid scan delay line operating at 20 Hz, with the pump amplitude 
modulated by a photo-elastic modulator (PEM) operating at 100 kHz. The reflected probe 
beam was detected by a silicon photodiode, which was output to a lock-in amplifier referenced 
to the PEM modulation. 

Measurements were performed on different cleaves of the same sample, comprising four 
data sets. The data corresponding to a Fermi surface angle of = 27° in Figs. 2, 3, and 4 are 
from the first cleave. The data corresponding to = 32°, = 38°, and = 45° (cut 5) in Figs. 
2, 3, and 4 are from the second cleave. The data corresponding to = 31° and = 45° (cut 6) 
in Figs. 1, 3, and 4 are from the third cleave. The time-resolved reflectivity data are from the 
fourth cleave. 
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ARPES Detector Non-Linearity Correction 



A common feature in even the best photoemission detectors is that the detector responds in a 
non-linear fashion to electron events (30-33). The effect is more pronounced with laser ARPES 
than synchrotron light because of the smaller background and higher electron flux. A typical 
curve showing the relationship between electron events and measured counts is displayed below 
in Fig. SI. 
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Figure S 1 : Calibration method for detector non-linearities. (A) Typical raw ARPES image, used 
for calibrating detector non-linearity effects. Images were taken at several different photocur- 
rents, and the detector response was measured in a high intensity region (Region I), moderate 
intensity region (Region II) and low intensity region (Region III). (B) Photocurrent vs. mea- 
sured intensities for the 3 regions shown in (A). Responses from the three regions have been 
vertically scaled to match a single curve. 

As can be seen from the figure, the response is linear for the lowest count rates (< 0.4 
counts/frame/pixel) but deviates from this at moderate count rates. In order to adjust for these 
effects, we measured the sample using several different laser intensities before proceeding to the 
pump-probe experiment. Using a picoammeter that is attached to the sample, we were able to 
associate each laser intensity with a photocurrent. Photocurrent and laser power are proportional 
to each other (34). By plotting photocurrent against the detector response in counts/frame/pixel, 
we then obtained the curve in Fig. SI, and we applied this curve to the pump-probe data to 
remove the influence of non-linearity effects in the final analysis. We note that the primary 
focus of the paper is the above -Ep spectral weight, where the count rate is low and thus non- 
linearity corrections are small. 
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Superstructure 



All ARPES data were taken using cuts through the Fermi surface along the T — Y direction in k- 
space in order to avoid complications with umklapp replica bands. The (y/2 x \f2)RA5° shadow 
band is beyond the fc-space windows examined and has no influence on the data. While higher 
order umklapp replicas of the main band and shadow band can be observed when a sufficient 
volume of statistics is accumulated, such replica bands are orders of magnitude dimmer than 
the main band at this photon energy, and are not observed in these data sets. 

Bilayer splitting in Bi2Sr2CaCu20 8+ g is visible at a photon energy of 5.9 eV, as seen in Fig. 
1 . In the off-nodal cut such splitting is visible in the form of the horizontal steak of intensity 
below the Fermi level for k > kp. Bilayer bands are also present along the nodal direction, 
although at the node they effectively lie on top of one another. 



Fluence Determination 

Assuming a Gaussian intensity profile, the total average fluence is given by 

Fave = (Geometrical factor) x 0.88 E pu 

FWHM^ + FWHM^' 

where E pu is the pump pulse energy, and FWHM pM and FWHM pr are the pump and probe full- 
width half-maximum intensity values. In the experiments here reported, the pump and probe 
spot size were measured in two different ways: a) by using a flip mirror to divert the beam away 
from the chamber, toward a CCD beam profiler; and b) by using a pinhole mounted directly on 
the sample manipulator. The two methods gave same value within 10% of the beam diameter. 
The geometrical factor arises from the fact that both pump and probe beams are incident on the 
sample from an oblique angle. When the beam approaches a sample at an angle 9 from normal 
incidence, it will spread into a spot that is larger by a factor of 1 / sin 9 . The largest correction 
resulting from this effect and applied to our data corresponds to a 20% adjustment downward 
in the fluence for the cut at § = 27°. 

In order to verify the decay rate momentum dependence and definitively separate it from 
potential systematic geometrical effects, two different geometrical configurations were used. 
Data corresponding to cuts 1, 3, 4, and 5 in Fig. 3 were obtained by rotating the sample flip 
stage, which has the benefit of giving sharp, relatively symmetric cuts across the Fermi surface, 
but results in a spot size that changes with Fermi surface angle for a given pump power. Data 
corresponding to cuts 2 and 6 were obtained by rotating the sample azimuth, where pump and 
probe spot sizes are identical between nodal and off-nodal cuts. The two types of measurements 
were in good agreement. 



Error Bars 

Special care was taken to determine the accuracy of the fluence measurements. The factors 
contributing to fluence error include variations in the power meter reading, losses from the 
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viewport windows, pump-probe misalignment, pump spot size, and probe spot size. Of these, 
pump-probe misalignment is the most important cause of statistical uncertainty between one 
fluence measurement and the next in the ultrafast ARPES experiment. A misalignment of 8 = 
0.2 x FWHMp H for a pump-probe ratio of 2:1 can result in an average incident fluence that 
differs from the 8 = value by 8%. We determined the combined fluence error according to 
the formula, 

(°F\ 2 = (°r\ 2 003 2 ( AS V 4FWH M^ g ,i+4FWHM^g p 2 ,- 
^ UJ- JTC. \FWHM P J (FWHMj l( + FWHMy 2 ' 

power meter ,. s v ' 

F alignment , . 

° spot sizes 

where F is the fluence, P is the laser power, A is a proportionality factor, 8 is the pump-probe 
spacial misalignment, and FWHM pu and FWHM pr are the respective pump and probe full- 
width half-maximum intensity values as in Eq. si. The variables o refer to the standard error 
of their respective subscripts. 

For the ultrafast ARPES measurements, we estimate 10% fractional power meter uncer- 
tainty, 3% fractional uncertainty due to window losses, 10% fractional uncertainty from pump- 
probe alignment, 10% uncertainty in the respective pump and probe FWHM values, and a typ- 
ical pump-probe ratio of 2:1, which leads to a combined fractional fluence uncertainty of 22%. 
For the ultrafast reflectivity measurements, the pump and probe had the same spot size, about 54 
jUm, and alignment was easier due to a smaller cryostat, comparable pump and probe spot sizes, 
and the fact that both beams were of the same frequency. We estimate 10% fractional power 
meter uncertainty, 3% fractional uncertainty due to window losses, 5% fractional uncertainty 
from pump-probe alignment, 10% uncertainty in the respective pump and probe FWHM values, 
and a typical pump-probe ratio of 1:1, which leads to a combined fractional fluence uncertainty 
of 18%. 

Horizontal error bars in Fig. 4 are based on potential skewing of the photoelectron exit 
angle from stray electric fields inside the chamber, based on a comparison between a tight- 
binding model for Bi2Sr2CaCu20 8+ g (35) and fc-space maps that were acquired at E = Ep at 
low temperature. An analyzer work function 4v = 4.35 eV was found to produce an adequate 
match between the data and the expected main band. We estimate the error in Fermi surface 
angle to be on the order of « 2° at <j) = 27°, and not more than Cty 1.3° in the nodal 
direction, which can be aligned more accurately because of the inherent symmetry of the band 
structure. An example map of the low temperature Fermi surface is displayed in Fig. S2. 

Vertical error bars in Fig. 3 are based on fit standard deviations. Vertical error bars in Fig. 
4 are based on the quadrature sum of the fit standard deviation from Fig. 3, and 22% of the fit 
value, which incorporates the uncertainty in fluence. 

Sample Heating 

An important systematic effect to rule out is the possibility that the fluence dependence is not 
indicative of the incident pulse, but is rather the manifestation of a slight rise in equilibrium 
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k x (7i/a) 



Figure S2: Low temperature constant energy map of Bi2Sr2CaCu20 8+( 5 at E = Ep. The ex- 
pected main band, in red, is from a tight binding model (35). 
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temperature caused by time-averaged heating due to the pump pulse. We tested for this by 
reducing the repetition rate of the laser from 543 kHz to 181 kHz while maintaining a fluence 
of 15 /iJ/cm 2 . The residual heating was thus equivalent to that caused by a fluence of 5 /iJ/cm 2 
at 543 kHz. Fig. S3 shows a comparison of decay rates as measured using the two settings. 
There is no discernible difference in the decay curves, which verifies that heating cannot be 
responsible for the observed fluence dependence. 
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Figure S3: Dependence of spectral decay on laser repetition rate. Red circles and orange dia- 
monds correspond to 543 kHz, which is the repetition ratio used for all data discussed in the 
main text. Blue and cyan triangles correspond to 181 kHz. 



Superconducting Gap Measurements 

The gap values in Fig. 2 were extracted according to the methodology of Norman et al. (72), 
which is a commonly accepted procedure in the field. We fit symmetrized EDCs at the Fermi 
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momentum to the convolution of a Gaussian (FWHM = 20 meV) and the function 



A(a>) = Const, x = . (s3) 

{(0-A 2 k /(o) 2 + ^ 

This amounts to a self energy ansatz of E(o) = A|/co — iT, where a phenomenological constant 
— iT has been added to the BCS self energy correction in order to account for finite lifetime 
effects. 

A graph of Eq. s3 is shown in Fig. S4A along with the resolution-convolved form that 
was used to fit the data. The relevant parameters are A^ = 19 meV and T = 10 meV, which 




Energy (eV) Energy (eV) Energy (eV) 



Figure S4: (A) Graph of Eq. s3 and its convolution with a Gaussian resolution function R(co) of 
FWHM = 20 meV. (B) Symmetrized EDCs and EDC fits at k F for a momentum cut at = 32°. 
(C) Symmetrized EDCs and EDC fits at kp for a momentum cut at = 27°. 

produce a curve that adequately matches the top symmetrized EDC in Fig. 2A of the main text. 
Because Eq. s3 has asymmetric quasiparticle peaks, the fitted gap size is slightly smaller than 
the maximum values of the resolution-convolved peaks would suggest. In Fig. S4B-C we show 
a magnification of the symmetrized EDCs from Fig. 2 of the main text corresponding to the 
highest pump fluence values. To complement the gap values extracted from the fitting analysis, 
we provide a comparison to gap values extracted from explicit measurements of peak locations 
in the symmetrized EDCs. At = 32° the equilibrium EDC peak at t = —0.3 ps is 18.6 meV. 
At t = 0.6 ps, this value is 11.6 meV, a reduction of 38%. At = 27° the equilibrium EDC peak 
shrinks from 27.9 meV at equilibrium (t = —0.3 ps) to 25.6 meV at t = 0.6 ps, an 8% reduction. 
These changes are in qualitative agreement with the change in gap size as measured using the 
fit, which is reduced by 42% at = 32°, and 19% at ^ = 27°. In both cases the gap at = 32° 
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is more strongly affected by the pump pulse than the gap at (j) = 27°. The fit yields a larger 
change in gap size than the peak-to-peak analysis, which may indicate that the finite energy 
resolution and data noise impose limits on the effectiveness of the peak-to-peak analysis, or it 
may reflect limitations of the fitting algorithm. Fig. 2 in the main text displays fitted gap values 
rather than results of a peak-to-peak analysis because the fit makes use of more data points than 
the peak-to-peak analysis, because it adequately captures asymmetries in the EDC lineshapes 
that would not be accounted for in simpler fitting algorithms, and because it rests on a firmer 
theoretical foundation than the simpler peak-based extraction method. 

Further comparison between the extracted gap and a simpler peak-to-peak analysis is shown 
in Fig. S5. The lighter triangles and circles correspond to the fitted gap parameter from Eq. s3. 
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Figure S5: (A) Normalized change in the equilibrium gap size for a cut at (j) = 32°. Cyan 
triangles correspond to the fitted gap parameter from Eq. s3. Blue inverted triangles correspond 
to 1/2 of the peak-to-peak distance of the fit curves in Fig. S4B. (B) Normalized change in 
the equilibrium gap size for a cut at = 27°. Light purple circles correspond to the fitted gap 
parameter from Eq. s3. Dark purple diamonds correspond to 1/2 of the peak-to-peak distance 
of the fit curves in Fig. S4C. 

The darker triangles and diamonds correspond to 1/2 peak-to-peak distances from the fit curves 
in Fig. S4B-C. As is the case in comparison to the raw data above, the fit parameter yields a 
larger change in gap size than the peak-to-peak analysis, but the cut at = 32° is more strongly 
affected by photoexcitation than the cut at = 27° under any analysis scheme. 
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Recombination Model 



Comparison of Reflectivity and ARPES Decay Rates. While the Rothwarf-Taylor model 
is a widely accepted framework for describing relaxation processes in the cuprates, there are 
competing views as to whether or not cuprate dynamics lie within the boson bottleneck limit of 
these equations. Under bottleneck conditions, hot bosons regenerate quasiparticle pairs much 
faster than they anharmonically decay or diffuse into the bulk of the sample. In the absence of a 
bottleneck, the opposite is true. Nevertheless, in the low temperature and low fluence limit both 
viewpoints reduce to the same functional form, which is the following at zero temperature, 

n ex =-Rn 2 ex n e - x 1 (t) = n e - x 1 (0)+Rt, (s4) 

or finite temperature, 

h ex = -R (n 2 ex + 2n ex n T ) n~ x \t) = ( + -*—) e 2n ? m - (s5) 

\2n T n ex {0) J 2n T 

Here n ex (t) is the population of photoexcited quasiparticles, n T is the population of thermal 
quasiparticles, and R is the intrinsic Cooper pair recombination rate in the bottleneck-free limit, 
or a renormalized Cooper pair recombination rate in the bottleneck limit. A more complete 
explanation of how these equations follow from the full Rothwarf-Taylor model is described 
elsewhere (17, 23). 

In Fig. 3 of the main text, ARPES and reflectivity quasiparticle response curves are com- 
pared at short times by using a fit to an exponential function f(t) = AIoe^ r °^^ to ^&(t — to). The 
short-time limits of Eqs. s4 and s5 are respectively 



n ex (t) 
n ex (0) 

and 

flprit) 



1-Rn ex (0)t + O(t 2 ) (s6) 



= l-R[n ex (0)+2n T ]t + O(t 2 ). (s7) 



n ex (0) 



In the case of ARPES, 6 eV photons induce photoelectrons with a mean free path of 5 
nm. This is well below the penetration depth of the 1.48 eV pump beam (100 nm), which 
means that only the surface of the material is probed. Furthermore, the initial photoexcited 
population n ex (0) vastly outweighs the thermal population nj in the majority of cases. To a 
good approximation, then, the decay rate fits displayed in Fig. 3 correspond to the product 
Rn ex (z = 0,f = 0). 

In contrast, the optical reflectivity probe penetrates the same distance into the sample as the 
pump, and therefore necessitates a more intricate analysis of spacial variation in n ex {z,t). The 
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problem has been addressed by Gedik et al. (17), who propose that the optical response should 
be modeled by inserting the solution of Eq. s4 into a weighted integral as follows: 



AR(t) = 



2AR(0)a 
2AR(0) 



POO 

/ dze- az n ex (z,t), 
Jo 

, ln(l + 7oO" 



n ex (z,t) 



n ex (0,0)e- az 
l+Rn ex (Q : Q)e- a H 



Jot 



(s8) 



(s9) 



The constant Jo = Rn ex (0, 0) corresponds to the initial decay rate of quasiparticles at the surface. 
Finite temperature effects are ignored because nj <C n ex (0) (as noted above). In the short-time 
limit, the decay rate of this revised equation differs by a factor of 2/3 from its counterpart at the 
sample surface: 

AR(t) 2 . ?N 

This is the most important difference between the physical interpretation of ARPES and optical 
decay rates in the short-time limit. 



Momentum-Dependent Recombination. In order to allow for momentum-dependent re- 
combination processes, Eq. s5 may be recast as follows: 

n k = — n k J Rkk' n k'd 2 k' — n k J R kk in k i T d 2 k' — n kT J R kk in k id 2 k' (sll) 

where n k is now the population of photoexcited quasiparticles at a given momentum k, R kk > is 
a recombination function related to the formation of Cooper pairs, and n kT is the momentum- 
dependent population of thermal quasiparticles. The first term on the right corresponds to pho- 
toexcited quasiparticles recombining with other photoexcited quasiparticles. The second term 
corresponds to photoexcited quasiparticles at k recombining with thermal quasiparticles in other 
parts of the Brillouin zone. The final term corresponds to photoexcited quasiparticles in other 
parts of the Brillouin zone recombining with thermal quasiparticles at k. The integral in the 
second term is a constant, so the term does not contribute to any density or fluence dependence. 
The third term vanishes for all momenta except the nodal direction, and along the node it is 
small relative to the first term on the right for all but the lowest excitation densities. At short 
time, then, momentum-dependent recombination dynamics are dominated by the first term on 
the right. 

In general, n k is a function of the fluence F, and can thus be written as n k (F). Because n k 
is the photoexcited quasiparticle population n k (0) = 0. As a result, for small fluence n k (F) can 
be expanded to linear order in F as n k ps 0C k >F. Making this substitution to the integrand in the 
first term results in 

n k ps — n k \ R kk m k id 2 k' ps — n k \ R kk i a k iFd 2 k' (s!2) 
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After pulling F outside of the integral and taking the derivative, the approximation finally 
duces to 

Yko ~ -F j R kk > a k > d 2 k' => -|? « j R kk , a k > d 2 k', (s 
where y k o = %o/%o is the momentum-dependent initial decay rate plotted in Fig. 3. 
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